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A complex quantum system can be constructed by coupling simple 
quantum elements to one another. For example, trapped-ionll^^l or su¬ 
perconduct ing®l quantum bits may be coupled by Coulomb interactions, 
mediated by the exchange of virtual photons. Alternatively quantum ob¬ 
jects can be coupled by the exchange of real photons, particularly when 
driven within resonators that amplify interactions with a single electro¬ 
magnetic mode®. However, in such an open system, the capacity of a cou¬ 
pling channel to convey quantum information or generate entanglement 
may be compromised®. Here, we realize phase-coherent interactions be¬ 
tween two spatially separated, near-ground-state mechanical oscillators 
within a driven optical cavity. We observe also the noise imparted by the 
optical coupling, which results in correlated mechanical fluctuations of 
the two oscillators. Achieving the quantum backaction dominated regime 
opens the door to numerous applications of cavity optomechanics with 
a complex mechanical system®®^. Our results thereby illustrate the po¬ 
tential, and also the challenge, of coupling quantum objects with light. 
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Cavity optomechanical systems comprised of a single mechanical oscillator inter¬ 
acting with a single electromagnetic cavity modeff^ serve useful quantum-mechanical 
functions, such as generating squeezed light^^^^®, detecting forces with quantum- 
limited sensitivity® or through backaction-evading measurement®, and both en¬ 
tangling and amplifying mechanical and optical mode^®. Systems containing sev¬ 
eral mechanical elements offer additional capabilities. In the quantum regime, these 
systems may enable two-mode backaction-evading measurements^^, creation of non- 
classical stated, fundamental tests of quantum mechanics®®, correlations at the 
quantum level with applications in high-sensitivity measurements®, and quantum 
information science®. Realizing these proposed functions requires multiple-element 
cavity optomechanical systems in which quantum mechanical optical force fluctua¬ 
tions dominate over thermal and technical ones. 



Figure 1: Schematic setup, a Two addressable oscillators are created by trapping 
ensembles of cold atoms in two sites of a superlattice potential in a cavity. The 
center-of-mass motion of each cloud of atoms is linearly coupled to the pump light, 
creating a system of two optomechanical oscillators, b Illustration of the realized 
system of coupled oscillators. Scattered pump photons induce coupling between 
the elements, and at the same time leave the cavity leading to backaction noise, c 
In a quantized picture, cavity photons act as the force-mediating particle for the 
interaction. Pump light is Stokes scattered off one element, generating a cavity 
photon that is absorbed through anti-Stokes scattering by the second element (and 
vice versa). Illustrations are not to scale. 
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An important new feature in these multi-mechanical systems is photon-mediated 
forces between mechanical elements. Consider a driven cavity containing two me¬ 
chanical elements with linear optomechanical coupling (Fig.[^). Each element ex¬ 
periences radiation pressure proportional to the number of intracavity photons. 
The displacement of one element changes the cavity resonance frequency, causing a 
change in the intracavity photon number, and thereby modifying the force on the 
second element. In this manner, an effective optical spring is established between 
the mechanical elements (Fig.[^). A quantized picture clarifies the role of cavity 
photons as the force-mediating particle for this interaction: pump light is Stokes 
scattered off one element, generating a cavity photon that is absorbed through 
anti-Stokes scattering by the second element, and vice versa® (Fig.[^). This cavity- 
mediated force has been shown to cause hybridizationllSHSl] and synchronization^ of, 
and incoherent transfer of energy betweenl^H^ mechanical modes. Cavity-mediated 
long-range interactions are also exhibited by the mechanical instability of an atomic 
gas within a transversely pumped optical cavity^^^^. 

However, in such an optically driven system, the cavity photons serve not only 
as force carriers but also as information carriers: The light emitted from the cavity 
permits measurements on the mechanical elements inside the cavity, leading to mea¬ 
surement backaction noise and limiting the capacity of cavity-mediated interactions 
to transfer quantum states between the coupled elements®. A present challenge is 
first to identify and then to ameliorate or to exploit this incoherent cavity-mediated 
force. Thermal and technical noise prevented the observation of such incoherent 
backaction in previous experiments. 

Here, we study a cavity optomechanical system containing two non-degenerate, 
spatially separated and well isolated mechanical oscillators, each comprised of a 
trapped ultracold gas of ®^Rb atoms. Applying cavity pump light that couples 
the mechanical elements, we observe light-induced forces through the redistribution 
of optomechanical cooperativity and, thereby, of radiation-pressure quantum noise 
between the different modes of motion. Additionally, we perform time-resolved 
measurements of the mechanical states of the two oscillators, with sensitivity near 
the standard quantum limit, to reveal the phase-coherent, oscillatory exchange of 
energy between the elements, and also the buildup of incoherent backaction noise 
that accompanies the cavity-mediated coupling. 

Following previous work^^, we trap around 900 atoms in each of two adjacent 
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wells of a superlattice potential, formed by two standing-wave optical fields that 
are resonant with TEMqo modes of a Fabry-Perot optical resonator and far detuned 
from atomic resonances. The deep superlattice potential prevents the exchange of 
atoms between the two ensembles. Optomechanical interactions are introduced by 
driving the cavity with light detuned by -42 GHz from the D2 atomic resonance. The 
driving frequency is close to another TEMqo cavity resonance, the frequency of which 
varies linearly with the center-of-mass displacement of each of the trapped atomic 
ensembles along the cavity axis. The linear optomechanical interaction strength 
is quantified by the optomechanical cooperativities of the two oscillators, given as 
Ci = Angf/{KTi), where gi is the single-photon/single-phonon coupling frequency 
for the oscillator, and k = 27 r x 1.82 MHz is the cavity half-linewidth. The 
transmission of this pump/probe light through the cavity is detected by a heterodyne 
optical receiver. At a typical average cavity photon number n = 2, Ci ^2 — 2 allowing 
for single-shot detection of the motion of each oscillator near the standard quantum 
limit^. Probing at this intensity reveals the resonance frequencies ({a;i,a; 2 } ~ 
27 r X {110,116.4} kHz) and linewidths (Pi ,2 ~ 27 r x 1.5 kHz) of the center-of-mass 
oscillations of the two ensembles (Supplementary Section SI). The superlattice 
potential allows us to tune the mechanical frequencies to be either non-degenerate 
(as here) or degenerate (see Methods). 

The motion of the two ensembles becomes coupled when the pump light is tuned 
away from the cavity resonance. In the unresolved sideband regime {k S> u;i^ 2 ); 
the optical force E) on oscillator i varies with the oscillator positions zi ^2 as Fi = 
— where the optical spring constants kij are given by 

, _ kiQigj 2nApc , , 

Zno,iZ}io,j + Ape 

Here, Ape is the difference between the pump light and cavity resonance frequencies, 
= \/h/2Mi(jJi quantifies zero-point motion and Mi is the total mass of the 
atoms in each lattice site (Supplementary Section S2). 

Owing to these optical forces, the mechanical normal modes, describing uncou¬ 
pled motion of each of the two oscillators in the absence of the optical spring, become 
increasingly coupled. As the optical spring gains strength, one normal mode evolves 
so that its position variable matches the sum of oscillator positions weighted by op¬ 
tomechanical coupling strengths. The optomechanical cooperativity of this “bright” 
normal mode becomes collectively enhanced. The remaining “dark” normal mode 
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Figure 2: Coupling of two shot noise driven oscillators. The cavity containing 
two distinct atomic ensembles is pumped with light at constant intracavity intensity 
(n ~ 5) and variable detuning Ape from cavity resonance, a The symmetrized 
sideband power spectrum of the transmitted pump light, calculated from a 5-ms- 
long measurement and averaged over about 250 experimental repetitions, is shown 
scaled to the shot noise power seen away from mechanical resonances. Spectra are 
offset by 3 shot-noise units from one another for clarity. Optical coupling of the 
oscillators results in a high-cooperativity normal mode at high frequency and a 
low-cooperativity mode at low frequency, b The coupling increases the frequency 
separation between the modes. Bars indicate standard fitting errors. Solid lines 
show theoretical predictions. 
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becomes increasingly decoupled from the cavity. The optical coupling also increases 
the frequency separation between the mechanical resonances. 

We observe both these effects, the renormalization and frequency separation of 
the mechanical modes, by monitoring the symmetrized mechanical sideband spec¬ 
trum of light transmitted from the cavity (Fig. [^. We examine the case Ape > 0, 
where the optical spring leads to mechanical stiffenin^^Sl^ constant intracavity 
photon number (n ~ 5), we vary the pump detuning, increasing the optical spring 
constant up to |A:i^ 2 | ~ 1 x 10“^^ N/m. We quantify the resulting coupling strength 
for our system by — 27r x 2 kHz, where for simplicity we approximate 

H by employing the average mass M of both oscillators (Supplementary Section S3). 
The steady-state spectra, with the illumination time much longer than the mechan¬ 
ical damping time, reveal the motion of the intracavity mechanical elements driven 
by zero-point fluctuations, thermal noise (at a level of around = 1-5 phonons), 
and radiation pressure quantum fluctuations. 

At zero detuning, we observe two uncoupled mechanical resonances with compa¬ 
rable line strength. At non-zero detuning, the frequency difference between the two 
mechanical resonances increases. The increase is relatively small because the optical 
coupling in our experiment is rather weak (H < a ;2 — wi); regenerative mechanical 
oscillations^ at higher optical power prevent us from increasing H further at steady 
state. The detuned pump also leads to a difference in the power of the two mechan¬ 
ical sidebands, with the higher frequency resonance becoming significantly stronger. 
The strength difference reflects not only that the mechanical modes are renormal¬ 
ized so that the cavity output is more sensitive to motion of the higher frequency 
normal mode, but also that this brighter mode has larger position fluctuations since 
it is more sensitive to radiation pressure fluctuations. A linearized model of cavity 
optomechanics incorporating both these compounding effects matches well to the 
data (see Supplementary Section S4). 

The cavity-mediated force should allow for the coherent exchange of energy be¬ 
tween the two oscillators. To examine this exchange, we perform two-dimensional 
spectroscopy, such as used to study vibrational coupling in molecules and solids®^, 
consisting of three steps. First, we excite a coherent state of motion (at a level of 
a few Zho — 0.8 nm) on just the low-frequency oscillator by resonantly modulating 
the intensity of one of the superlattice trapping beams. During this time, the cavity 
is probed with on-resonant light (h = 2) so that the mechanical oscillators are not 
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Figure 3: Phase-coherent energy exchange, a The optical cavity probe mea¬ 
sures a weighted sum of the displacement of the two mechanical elements. After 
one element is coherently driven and an optical link between the modes is estab¬ 
lished, the beat note in the mechanical ring-down signal shows coherent oscillations 
of both elements. Data (black) for each Tc, averaged over about 10'^ experimental 
repetitions, show reasonable agreement with zero-free-parameter numerical simula¬ 
tions (red), b Phase-space plots indicating the average relative state of the high- 
(blue) and low-frequency (green) oscillators just after applying the optical coupling 
for varying times Tc show two cycles of energy exchange. Numerical simulations 
take into account experimental uncertainties (red shaded area) and show reasonable 
agreement. Bars denote standard statistical errors. 


coupled. Second, we turn on the optical coupling between the oscillators, at the 
level of D ~ 27r X 4 kHz, by tuning the pump/probe light away from resonance 
(Ape = 27r X 1.4 MHz) and increasing its intensity (h ~ 8), holding the light at 
these settings for a variable time Tc- Third, we turn off the optical spring and record 
the phase modulation of resonant cavity probe light (h = 2) to characterize the 
mechanical state (Supplementary Section S5). 

Averaging the phase modulation signal over many experimental repetitions, we 
observe a time-domain beat note in the mechanical ring-down signal (Fig.j^) show¬ 
ing that the coherent excitation is now shared by both oscillators. Applying matched 
data filters to this time trace, as demonstrated in Ref. [30|for a one-oscillator system 
and detailed in the Supplementary Information (Section S6), we obtain measure¬ 
ments of both quadratures of motion of each mechanical element at the time when 
the optical spring is turned off (Fig. [^). We observe a cycle time for the exchange 
of energy between the oscillators of about 100 fj,s. The maximal transfer of energy 
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occurs after a half-cycle of oscillatory exchange, at the level of A^I{A\ + A^) ~ 0.2 
where Ai is the excitation amplitude of oscillator i. We compare these quadrature 
measurements to numerical simulations (Supplementary Section S7). The agree¬ 
ment with data is reasonable given the effects of systematic uncertainties in h and 
Ape, and fluctuations in M* and gi. 

Our experiment demonstrates that two mechanical elements with no direct con¬ 
tact can be coupled to one another controllably through cavity optomechanics. We 
identify two conditions required for this controlled spring force to induce fully co¬ 
herent quantum evolution of the coupled mechanical system. First, one requires the 
coupling rate to be larger than the thermal mechanical decoherence rate (r'th + l)rM. 
This condition is marginally satisfied in our work. Second, one requires coherent 
evolution of the optical field acting as the force carrier. In our system, the non-zero 
cavity linewidth implies that the field leaks from the cavity, providing a means to ob¬ 
serve the state of the optical spring. The force noise resulting from this observation 
limits the coherence of the mechanical system. 

We observe the consequence of this force noise by applying the aforementioned 
matched filters to each individual repetition of the three-step experimental sequence 
described above (Fig. |^. Measurements before the oscillators are driven and cou¬ 
pled characterize their initial state. The distribution of quadrature measurements 
is centered about zero. The excess variance in the quadrature measurements above 
the optical shot-noise level reflects the incoherent initial excitation of both oscilla¬ 
tors at the level of 2.5 phonons. After the oscillators are coherently excited and 
optically coupled, the quadrature measurement distributions are offset from the ori¬ 
gin, reflecting the coherent motion of both elements. Additionally, the variance of 
these measurements increases, showing that backaction noise during the application 
of the optical spring incoherently excites both oscillators. The incoherent excitation 
strength is different for the two oscillators, with the high-frequency oscillator gaining 
energy faster than the low-frequency oscillator. This difference occurs because the 
bright coupled mechanical mode is more strongly disturbed by backaction, and, in 
the weak-spring regime explored in this work, this bright mode corresponds to larger 
displacement of the high-frequency oscillator than of the low-frequency oscillator. 
Our observations are in good agreement with expectations from zero-free-parameter 
theory (see Supplementary Section S7). 

Interestingly, the backaction noise results in correlated fluctuations of the two 
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Figure 4: Buildup of correlated backaction noise during optical coupling. 

a Distributions of about 10^ quadrature measurements on the two oscillators under 
three conditions: no coherent excitation or coupling (left) or coherent excitation and 
coupling of duration Tc (center and right). Dashed lines indicate l-cr variance ellipses 
for signal (red) and optical shot noise (black); the excess signal variance quantifies 
the incoherent excitation. Optical coupling increases the incoherent excitation level, 
plotted in b as the number of added phonons for the high- (a) and low-frequency 
(t) oscillator. The measured data show agreement to zero-free-parameter theory 
calculated for each Tc (open symbols); see Supplementary Information. The coupling 
leads to correlations in the displacements of the two oscillators, quantified by the 
covariance r (♦) between their quadratures. Correlations show the expected trend, 
slightly reduced below zero-free-parameter expectations (open symbols). The mea¬ 
sured 1-a correlations for uncoupled oscillators are indicated by grey shaded range. 
Solid lines indicate theory for average experimental parameters. Bars show typi¬ 
cal l-cr standard statistical errors; a plot with all error bars can be found in the 
Supplementary Information. 
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oscillators. These correlations, quantified by examining the covariance between 
motional quadratures of the two oscillators, grow in strength during the application 
of the optical spring. We clearly observe the build-up of correlations that is larger 
than the measured level of uncertainty for uncoupled oscillators, and is in reasonable 
agreement with zero-free-parameter theory (Fig. |^. 

Fluctuations of the optical spring force are found to add one phonon of inco¬ 
herent excitation within about 30 fis, a time shorter than the coherent coupling 
time (0/27r)“^ ~ 250 fis. This comparison implies that for the conditions of our 
experiment - pump light detuned about one half-linewidth from the cavity reso¬ 
nance - the force linking the two mechanical oscillators cannot be used to exchange 
non-classical states of motion, such as phonon Fock states, between the oscillators. 
However, this limitation is not fundamental, and may be mitigated by driving the 
cavity further from its resonance frequencj^, or through post-selection or feedback 
using light exiting the cavity. 
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Methods 


Preparation of optomechanical oscillators. Sample preparation is similar to 
the procedure described in Ref. |25l After laser cooling and subsequent evaporative 
cooling in a magnetic chip trap, a sample of ®^Rb atoms is magnetically transported 
into the cavity. The atoms are then first transferred into a (single color) optical 
lattice, created by exciting a TEMqo mode of the cavity at a wavelengths of Aa = 843 
nm, with degenerate oscillation frequencies in all lattice sites. To create separately 
addressable oscillators, we excite a second TEMqo mode at Ab = 862 nm. The 
resulting superlattice (illustrated in Article Fig. 1) induces an adjustable frequency 
splitting between neighboring sites, which we tune in this work to be around 6.4 
kHz. We perform a final stage of evaporative cooling from the optical superlattice 
trap. The thermal phonon occupation in the superlattice (with the probe light 
off) is calibrated by employing time-of-flight measurements of the atomic cloud to 

~ 1.5 phonons in the relevant axial direction of the lattice. 

The atom population of the superlattice wells is controlled by carefully position¬ 
ing the magnetic trap before ramping up the optical lattice. The number of occupied 
wells is adjusted by the spatial extent of the cloud, tuned by rf-evaporation. For 
this work, we typically load each of two sites of the optical lattice with about 900 
atoms. The total number of atoms A^a is derived from the average shift of the cavity 
resonance frequency Awc = Wa5o/(2Aca) due to the presence of atoms, where qq is 
the single-atom cavity QED coupling and Aca the detuning of the probe light from 
atomic resonance. 

Optomechanical interaction. Information about the motion of both oscillators 
is extracted by optomechanical interactions with the probe beam, detuned by 
Aca — ~42 GHz to the D2 transition [kp = 27r/(780 nm)). The optomechanical 
interaction depends on the gradient of the probe light at the position z of the 
atoms, quantihed by sin(2fcpz). Choosing the two optical lattice sites closest to the 
maximum probe light gradient, we create two oscillators with comparable (nearly 
maximal) linear optomechanical interaction strength. The non-zero curvature of 
the probe light causes residual quadratic coupling®. We measure this additional 
contribution to the trapping potential by the probe independently to be about 
[wqp ~ — 27r X 40 (a;q^2 — +2vr x 208)] Hz/photon for the low (high) frequency 
oscillator. The fact that one oscillator experiences trapping whereas the other 
anti-trapping by the pump/probe field confirms the loading position around the 
point with maximum linear optomechanical coupling. 
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Supplementary Information 


SI: A gas of cold atoms as optomechanical oscillator 

With the preparation described in the Methods, we create a system where the 
motion of the trapped atoms is linearly coupled to the intracavity pump light. The 
cavity senses the collective center-of-mass mode of this motion along the cavity axis, 
whereas all other modes average out and are not measured. The collective center- 
of-mass motion of the cold gas can be thought of as equivalent to, for example, the 
motion of a SiN membrane and other realizations of optomechanics. This analogy 
establishes the description of our system as an optomechanical oscillator^®. We 
note that the equivalence of different realizations of optomechanical systems is 
reflected in a common theoretical description (see Ref. S2, for example). 

Extraction of oscillator properties. Fits of the mechanical sidebands observed 
in the phase quadrature of the optical heterodyne detection when probed on cavity 
resonance are employed to extract the mechanical damping rates Tj, resonance fre¬ 
quencies uji and the peak height of the incoherent response P^. The optomechanical 
cooperativity Ci is quantified for each oscillator by 

a = + 1/2) + + + (1) 

where 5 sn is the total shot-noise power spectral density and e 0.05 is the 
efficiency of the heterodyne detection^lSl. The single-photon/single-phonon coupling 
strength gi for each oscillator is then determined by Ci = Angf/{ kT i), with the 
mean intracavity photon number n. In this work, we typically prepare samples 
with <7* ~ 27 r X 24 kHz. 


S2: Optical forces and induced optical spring 

The optomechanical interaction of the cavity light with the oscillators can be de¬ 
scribed by an optical dipole force acting onto oscillator i which is linear in intracavity 
photon number n: 


Fi = - 


figi - 
-n, 


ZHO,i 


( 2 ) 


with ZHO,i the harmonic oscillator length. At the same time, a displacement Zj of 
an oscillator from its equilibrium position induces a change in the cavity resonance 
frequency by an amount gjZjjZnoj- In & cavity driven with a constant frequency 
pump, this change in resonance frequency translates to a change in intracavity pho- 
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ton number n, which we approximate for small cavity resonance frequency changes 
by 


5n ~ 


2Apc gj 

— Zj- 

+ ^pc ^HO,i 


( 3 ) 


With this, the force on oscillator i is given by 


« = E 

j 


^9j9i 2?T.Apc ^ 


^i,j Zj , 


( 4 ) 


where we have dropped constant terms and introduced the spring constant of Article 
Eqn. 1: 


^9i9j 2fiApQ 

Zno,iZno,j 


S3: Simplified classical model 

While the individually measured parameters are used to analyse the coupled os¬ 
cillator system for each experimental data set corresponding to a specific Tc, it is 
illustrative to consider a simplified system which also helps to quantify the coupling 
strength 0,. For equal spring constants k = kij, but different natural frequencies 
LOi and masses Mj, the classical equations of motion for our system of two coupled 
oscillators are given by 

_[ Ml ^ Mi{uj + 5/2)^+ k k \ f 

V 0 M2 ) \h) V k M2{Cd-5/2y + k J yz2j ' ^ ^ 

with the mean resonance frequency to = (wi + a; 2)/2 and the natural frequency 
splitting 6 = L 02 — uji. For equal masses M = Mi^ 2 , this further simplifies to 

{Q + 5/2)2 ^ X A 

k/M {oj - 6/2)^ + k/M J \z2) ' ^ ^ 

Normal mode frequencies of the coupled oscillators are given by the eigenvalues of 
the matrix in Supplementary Fqn. Q: 

Uj\ = 52/4 + k/M + w2 ± ^/{k/MYT^. (7) 

The system is in the weak (strong) coupling regime, if the splitting due to cou¬ 
pling is small (large) compared to |5|. Corresponding to Supplementary Fqn. Q, 
we quantify the strength of the coupling by 12 = In this work, we typically 







operate at O ~ 27r x (2 — 4) kHz < <5 = 27r x 6.4 kHz, i.e. in the weak coupling 
regime. 

The average oscillator properties for the entire data set used in this work differ 
by less than 15 %, which justifies the use of the simplified coupling strength H in 
the discussion of general features of our system. Calculations performed to predict 
the experimental observations take variations of gi, Mi and Tj into account (see 
below). 


S4: Linear amplifier model extension to two oscillators 

An optomechanical system can be modelled as a linear amplifier. For a single 
oscillator, this treatment was discussed in detail in Ref. S2. In order to model 
the steady state spectra for two oscillators optically coupled to the same cavity field 
(as presented in Article Fig. 2), we extend this model by allowing two oscillators to 
interact with the cavity field. In the framework of this treatment, the evolution of 
the field operators in frequency space can be expressed by 
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where Fa is the response matrix for the optical cavity, Fbj the matrix corre¬ 
sponding to the mechanical response of oscillator j and Tj the matrix capturing 
the coupling of oscillator j to the light field of the cavity. Photonic and phononic 
inputs (including noise) are symbolized by ain+, ain- and respectively. 


S5: Coherent excitation and optical coupling 

Resonant coherent drive. We study the coherent drive of oscillators at Ape ~ 
0, h ~ 2 and therefore ki 2 = 0, i.e. without coupling. 

First, we intensity modulate one of the superlattice beams with a rectangular 
pulse of only 3 cycles duration at Wexcite = 27r x 110 kHz. Because of the correspond¬ 
ing spectral width of about 37 kHz, such a pulse is resonant with both oscillators. 
In Supplementary Fig.[^ the resulting ringdown signal is presented. The observed 
signal shows the beat note of the two oscillators. For the time domain plots, the 
data are filtered by a second order butterworth bandpass filter to remove technical 
noise and shotnoise, (forward and backward, to preserve phase) with a bandwidth 
of 40 kHz, centered around the oscillators’ resonance frequencies. 
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Addressing only a single oscillator in a short time requires tight control over the 
resonance frequencies. The intensities of both lattice beams are detected after trans¬ 
mission through the cavity and actively stabilized. However, drifts in the detection 
efficiency cause the center frequencies of both oscillators to slowly vary over time. 
We therefore directly stabilize the low frequency oscillator to wi = 27r x 110 kHz, at 
Ape ~ 0, h ~ 2. For this, we extract the deviation of oji from a sideband fit averaged 
over a few experimental repetitions. A digital proportional-integral controller feeds 
back to the intensity of one of the superlattice beams, thereby stabilizing the center 
frequency of the low frequency oscillator (bandwidth on the order of a few minutes). 
The frequency drift of the second oscillator is observed to be limited to about ± 200 
Hz, since we only apply feedback to one beam of the superlattice trap. With the 
typical mechanical damping rate of ri^2 ~ 27r x 1.5 kHz, both oscillators are stabi¬ 
lized to within a fraction of their linewidth, allowing for stable spectral addressing. 
The oscillators’ frequencies are locked in this way for all experiments presented in 
this paper. 

For single oscillator addressing, we resonantly excite the low frequency oscillator 
with a pulse at Wexdte = 27r x 110 kHz and a duration of 50 cycles. We shape the ex¬ 
citation pulse with a truncated Blackman envelope to suppress spectral crosstalk® 
The small spectral width of such a pulse ensures the residual excitation of the high 
frequency oscillator with l <;2 ~ 27r x 116.4 kHz to be on the order of 10“^. Supple¬ 
mentary Fig.j^shows the resulting ringdown signal, dominated by a single oscillator. 

Coupling pulses. To turn the optical spring on and off, the probe/pump settings 
are varied gradually over 20 ^s to prevent spurious excitation of the mechanical 
elements from the increase in n. 

To facilitate data taking, we repeat excitation and coupling pulses in each single 
run of the experiment. After applying a pulse, the mechanical system is allowed 
to re-equilibrate under constant illumination (at n ~ 2 and Ape — 0) before the 
spectroscopy sequence is repeated after 2 ms. After 5 repetitions, the atoms are 
released from their trap and a new atomic gas is cooled and placed within the 
cavity. 

Choosing a sufficiently slow repetition rate allows to read out different parts of 
the timetrace: The first time-window covers the ringdown following the excitation 
and coupling. The second time-window is shifted by several ringdown times, so that 
any coherent excitation has damped out and each oscillator is exclusively driven 
by optical shot noise and in thermal equilibrium with the phonon bath. This time 
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window provides a record of the intrinsic properties (i.e. without coupling/drive) of 
each oscillator. The temperature of both oscillators in this timewindow corresponds 
to the initial incoherent excitation. We measure this initial phonon occupation 
for each Tc value and each oscillator by sideband asymmetryto be typically 2.5 
phonons. With an optomechanical cooperativity of Ci ^2 ~ 2, corresponding to 
about one phonon due to backaction®, and a thermal base occupation of about 1.5 
phonons (see above), we find good agreement with expectations. 

Coherent energy exchange rate. The effective frequency splitting between the 
two oscillators, determining the coherent energy exchange rate, is increased by two 
effects during the coupling. First, the coupling itself increases the natural frequency 
splitting of 6/2tt ~ 6.4 kHz by an amount of about 1.5 kHz. Additionally, the 
trapping potential of the pump light further adds about 1.5 kHz to the frequency 
difference due to the higher photon number during coupling (see Methods). We 
therefore expect a cycle time for coherent exchange of energy of around (9.4 kHz)“^. 
Our observations are in reasonable agreement with this prediction (Article Fig. 3). 

S6: Matched filter 

We employ a matched filter to extract quantitative results from timedomain signals. 
First, we extract Fi ^2 and oji ^2 from shotnoise-only driven oscillators for each Tc 
data set. We define two template functions ej{t) = The optimal filter 

coefficients aj are defined as the best choice to express a signal timetrace s{t) in 
terms of the template functions, i.e. 

s{t) = J2ajejit). 

3 

The inner product between a single timetrace s{t) and the template functions mea¬ 
sures their correlation and is given by 

ek{t)s{t)dt 

- couple 

where tend is the end of the ringdown timewindow. Plugging the expansion of s{t) 
into this expression gives 

{ek{t),s{t)) = 

3 

with Ekj = {ekit), ej{t)). The optimal filter coefficients can now easily be extracted 
by the multiplication of the vector {ek{t),s{t)) with the inverse of the matrix Ekj. 
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This procedure removes the effects of the template functions not being orthogonal. 
The aj are exactly the coefficients which minimize the accumulated deviation of 
the timetrace s{t) from the expansion Ylj constitute the best estimate 

of the state of each oscillator in phase-space at Tq- A direct comparison between 
the distribution of filter coefficients and timetraces is given in Supplementary Fig.[^ 
where we compare a time domain fit with amplitudes and phases as free parameters 
to the corresponding distributions of filter coefficients. The resulting phases and 
amplitudes agree with the matched filter results to within 10“^. Additionally, the 
filter coeofficients contain information about the noise of each oscillator through the 
phase-space area of its distribution. 

We further scale both real and imaginary parts of the filter output to harmonic 
oscillators units by multiplying with 1/y/eSsNCiTi (Ref. [T6l) . This procedure 
corrects for differences in Tj and Ci of the two oscillators and allows us to identify 
the filter coefficients with dimensionless harmonic oscillator displacements Zj. 

S7: Numerical simulations for pulsed experiments 

Using numerical integration, we calculate the expectation for the coherent energy 
exchange and the backaction noise due to coupling. 

The lab-frame master equation for the density operator p describing the cavity 
field and the mechanical modes including optical pump and dissipation is given by 

d _■ 2 

+ p]+ + ^+ ^)^{bj)p+ ^jl^th,j^{b^j)p ) 

i=i 

( 10 ) 


Tim = '^f^jb]bj ( 11 ) 

i=i 

Ho = biOcO^a -|- -|- ( 12 ) 

2 

T-(-om = ^ hgja^a{bj -h &J) (13) 

Ci6)p = dp6^-^(6^dp +po^o^, (i4) 


where we have introduced the cavity resonance frequency cuc, the optical pumping 
rate p and frequency ljl, the mechanical (optical) dissipation rates rj(K) as well as 
mechanical reservoir occupation numbers I'thj- 
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Figure 1: Resonant coherent excitation without coupling, and comparison 
of time domain fits and matched filter, a Time domain data (black) after 
exciting both oscillators (upper panel) and a single oscillator (lower panel). Blue 
shaded region indicates the timewindow in which the coherent drive is applied. 
Orange solid line shows the time domain fit in the ringdown timewindow. b Phase- 
space plots for the same data generated with the matched filter, in arbitrary units, 
after exciting both oscillators (upper panel) and a single oscillator (lower panel). 
Dashed red line shows lu-error ellipse generated from the data. Extracted parameter 
for each approach are given. Each data set contains about pulses. 
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The coupling of each oscillator to its independent, Markovian bath is captured 
by the third term of the master equation ( [lo) ). This term leads to dissipation and 
decoherence of oscillator j with a rate + towards a thermal state with t'thj 

phonons. Using this general formulation avoids assuming a particular dissipation or 
de-phasing mechanism or a preferred quadrature of motion. We chose this model 
because there are likely a number of microscopic relaxation and dephasing mecha¬ 
nisms and their details are not known. The good agreement of the experimentally 
obtained ringdown traces with the theoretical predictions supports our choice. 

Phonon annihilation and creation operators are related to the displacement op¬ 
erator by 


Zi 



(15) 


and the dimensionless, complex phasor Z used in the main paper is identified 
with -^{b) for each oscillator. Moving to a co-rotating frame for the mechan¬ 
ical modes, hi —>• linearizing the interaction and eliminating the cavity 

field yields the master equation for the reduced density operator p^. Introducing 

b = (6i, 62, ^1, ^2) 


where M, are the 4x4 matrices 
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For our system, we can safely discard terms oscillating at + LJ 2 ) and separate 
the remaining terms into unitary and non-unitary parts to arrive at 


dpr, 


dt 


^ [n,Pm\ + \ + ^^jP)) > 

*d=i ^ ^ 


with the effective Hamiltonian 

n = 6nib\bi + 5n2blb2 + + b\b 2 e^^‘^^-^^^*), 


where the frequency shifts and coupling strength are given by 

5Qi =Im {[0^]ii - [O 2 U) 

J =lm{[Of]i2 - [O^hi) ■ 


(24) 


(25) 


(26) 

(27) 


Simulation of motion in the time domain. We model the coherent motion of 
the coupled oscillators by the classical equations of motion, taking into account the 
specific parameters {gj,fi, u)j, Ape, t'thj, Fj, Mj} for each Tc data set and for each os¬ 
cillator determined as described above. We further take into account the additional 
frequency shift due to residual quadratic coupling, see above. To account for experi¬ 
mental imperfections in the created coupling pulses (slightly varying transients), we 
feed the averaged, measured waveforms for h and Ape for each Tc data set into the 
simulation. We obtain the time domain evolution of motion of both oscillators, as 
shown in Article Fig. 3, by numerical integration. 

To estimate the uncertainty of the predicted evolution, we vary the experimental 
parameters according to their uncertainty by randomly choosing values from a cor¬ 
responding Gaussian distribution. Repeating the simulation for a few hundred sets 
of random samples gives a phase-space area indicating the uncertainty, as illustrated 
in Article Fig. 3 as red shaded areas. 

We note that shot-to-shot fluctuations of the loading position and therefore 
of the atom number in each oscillator would decrease the coupling strength and 
therefore the observed coherent energy transfer. However, the single shot resolution 
of our experiment does not allow a quantitative estimation of this effect, and the 
agreement with expectation suggests this effect to be small. 

Model for dissipation. To calculate expectations for the backaction noise, we 
integrate the master equation using independently measured oscillator parameters 


X 



for each holding time Tc- The matrix h describes the non-unitary evolution of the 
mechanical degrees of freedom. It results from their coupling to independent heat- 
baths and the additional correlated noise as a consequence of the coupling mediated 
by the dissipative cavity field. Its explicit form is 
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and the parameters read in the unresolved sideband regime 
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The lines in Article Fig. 4 are results of numerically integrating the master equation 
Supplementary Eqn. (24) for the second order expectation values employing the 
experimental parameters averaged over all Tc and including the ramps to establish 
and turn off the optical spring. Open symbols in Article Fig. 4 are calculated by 
integrating the Supplementary Eqn. (24) numerically for independently measured 
experimental parameters h, ujj, Ape, t'thj, Fj}. 


Calculation of added phonons. The phase-space distribution measured by the 
matched filters is employed to determine the number of added phonons. Error 
ellipses give the phase-space areas Aj-ef for shotnoise driven oscillators and Aj-d after 
excitation and coupling, the contribution of shotnoise is accounted for by calculating 
the output of the filter off mechanical resonance (black error ellipses in Article 
Eig. 4). The number of added phonons (Article Fig. 4b, top panel) is then given by 
Az^i ,2 = , where (P_) is the power of the blue (red) sideband. 

We estimate the standard statistical errors by error propagation from uncertainties 
in sideband power and phase-space areas. We show typical uncertainties in Article 
Fig. 4 and plot the full results in Supp. Fig.[^. 

We test the statistical agreement between theory and experiment by calculating 
Tred ~ = 0.97 « 1, where N is the number of observations. 

The near unity result indicates a good statistical agreement between theory and 
experiment. The weak oscillatory behaviour of the added phonons predicted by 
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Figure 2: Added phonons and correlations plotted with all errors, a Article 
Fig. 4b, top panel, plotted with all l-ci statistical errors, b Article Fig. 4b, lower 
panel, plotted with all 1-a statistical errors. 


theory for average parameters cannot be resolved due to statistical uncertainties 
and to fluctuations of the experimental parameters between different Tc runs. 


Noise correlations. Noise correlations are quantified by 

(Re(ii)Re(i2)) 


’’Re = 


(Re(Zi)2)(Re(Z2)2) 


(31) 


and similarly for the orthogonal quadrature Im(Zj). Theoretical predictions are 
calculated by integrating Supplementary Eqn. ( [^ . For our system and conditions 
the theory predicts non-zero correlations r^e = Hm = We can thus treat the 
correlations between real and imaginary parts of the coefficients as independent 
measurements of the correlation r. To obtain correlations from the measured filter 
coefficients we account for non-ideal detection efficiency and correlations introduced 
due to the application of matched filters. The former effect is included by subtract¬ 
ing the variance of filter coefficients obtained far off mechanical resonances from the 


variances appearing in the denominator of Supplemental Eqn. (31). The correlations 


introduced due to the filters are obtained by subtracting correlations of uncoupled 
oscillators and shot-noise in the numerator. The thus corrected correlation coef¬ 
ficients can be directly compared to coefficients calculated by integrating Supple- 
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mental Eqn. ( |24| ) for independently measured experimental parameters for different 
coupling times (open symbols in Article Fig. 4b and Supp Fig. 2b). Statistical errors 
are obtained by performing a Fisher transformation on the statistical variance of 
the sample. To verify the validity of the statistical analysis, we also extract the co¬ 
efficients of the uncorrelated quantities (Re(Zj)Im(Zj)), which are scattered around 
zero with standard deviation given by the shaded region of Article Fig. 4b. 

The observed correlations, averaged over all coupling times, fexp ~ 0.152 are 
slightly reduced from zero-free-parameter theory expectations rtheo ~ 0.176. This 
deviation is probably caused by systematic uncertainties or shot-to-shot fluctuations 
leading to inhomogeneous averaging (see above). The oscillatory behaviour of the 
correlations predicted by theory for average parameters can not be resolved due to 
fluctuations of the experimental parameters between different Tc runs. 
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